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Abstract 

We determine the ground state energy of atoms and quantum dots whose number N of electrons 
is large. We show that the dominant terms of the energy are those given by a semiclassical Hartree- 
Fock theory. Correlation effects appear at the order N\nN for atoms and the order ./V for quantum 
dots. We compute them. The semiclassical Hartree-Fock theory creates oscillations in the ground 
state energy as a function of N. These oscillations reflect the dynamics of a classical particle 
moving in the presence of the Thomas-Fermi potential. The dynamics is regular for atoms and 
some dots, but we present the case of a dot where this dynamics is fully chaotic and we compute 
the oscillating part of the ground state energy in this case. 
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Since the introduction by Thomas [l| and Fermi [2j of the statistical theory of atoms, 
many efforts have been made to systematically improve on it. For neutral atoms with N 
electrons, an expansion of the ground state energy, whose first term is Thomas-Fermi (TF) 
theory, was determined up to the order N 5 ^ 3 . This expansion was based on Hartree-Fock 
(HF) theory. But after Lieb and Simon |4( proved that TF theory becomes exact when 
iV — > oo, part of the expansion was put on a rigorous basis [5|. A crucial step further, inside 
HF theory, was made by Schwinger and Englert who showed that there exist oscillating 
corrections to the ground state energy, of order iV 4//3 , and period iV 1 / 3 . They interpreted 
such terms as indicating shell effects. 

We consider quantum dots as artificial atoms made of N electrons, confined by a potential 
and constrained to move in two dimensions. The determination of the dependence of the 
ground state energy of quantum dots as a function of iV is of great interest. It was proven 
by Lieb, Solovej, and Yngvason [?| that in this case also, TF theory becomes exact when 
N — > oo. This suggests that like in the case of atoms, corrections to it should be determined. 
An important difference with atoms is the fact that the confining potential isn't necessarily 
rotationally symmetric, and can be changed. 

We have found a general formula for the ground state energy of both atoms and dots, 
appropriate to the case when iV -1 is small. It is based on the fact that in this case, we 
are confronted to a semiclassical mean field theory corrected by fluctuation effects. In this 
way, we have shown that HF theory is valid, up to a certain order in iV -1 where correlation 
effects become important. We have determined the dominant correlation correction, which 
in the case of atoms, is NlnN. 

We show that the oscillating contributions to the energy appearing in HF theory can be 
interpreted as coming from periodic orbits of a classical particle moving in the TF potential. 
In atoms, this motion is integrable. But in dots, the confining potential can be chosen so 
that the motion is fully chaotic, and we compute the ground state energy for a specific 
example of such a chaotic quantum dot. 

We can put non-relativistic atoms and quantum dots on the same footing, if we describe 
them, in appropriate units, by the d- dimensional Hamiltonian of N electrons 



For atoms (d — 3), Z is the number of protons and V(x) — — \x\ 1 . For dots (d — 2), Z will 
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be identified with N, a possible dependence of V(x) on N being neglected at this stage. 

For the purpose of considering systems where Z (and N) is large, we rescale the coordi- 
nates so that H = Z 2 ~ 2 / d H, where 

N N N 

tf^E^+Ew^E^, (2 ) 

j=i j=i j^k ■> 

the parameter e, playing the role of h, being e = Z~ x l d . In this form, the Coulomb interaction 
is of mean field type, and, when Z is large, the Hamiltonian has to be considered in the 
semiclassical limit. Our strategy to determine the ground state energy E of H is then the 
following: considering first e and Z as independent, we compute an asymptotic expansion in 
Z^ 1 of the pressure P of the Hamiltonian H, then take the zero temperature limit (/3 — ► oo), 
and finally look at the semiclassical limit e — ► 0. 

With this strategy in mind we first make a Hubbard-Stratonovich transformation on the 
Coulomb interaction term in the coherent state path integral representation of the partition 
function Q [8]. After an appropriate shift, the partition function can be written as the 
average value over an auxiliary Gaussian field 9 

Q={e Zm ), (3) 

S(8) being an action such that its saddle point, which dominates Q when Z is large, corre- 
sponds to the Hartree approximation. 

The final exact result for E can be put in the form 

E = E H + E X + E c;1 + E c . 2 . (4) 

The first term is the Hartree energy, given by 

E^iiN-%1 -^L po ( x )p (y)-2 r<leD{e). (5) 
2 J \x-y\ J 

D(e) = tr?9(e — h), with d the Heaviside step function, is the integrated density of states of 

the one-body Hamiltonian h = — e 2 A + W(x), with W(x) the self-consistent potential given 

by 

W(x) = V(x) + / -J^—My). (6) 

J \ x y\ 

Po(x) is the density po(x) = 2n(x\x)/Z , where n(x\y) = (jb|i?(// — h)\y) is the density 
matrix associated to h. Finally, the chemical potential fi is fixed by the constraint d^E = 0. 
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The second term is the exchange energy 



1 f dxdy o 
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This term also appears in HF theory, where it is the dominant exchange contribution when 
Z is large [9]. The remaining terms describe correlation effects, with 
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and 



E c . 2 = lim — \ ln(exp 
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C is the operator of kernel 



(7(03,^12/,^) = | — - — -5{r x -T y ), (10) 



\x-y 

< t x , < /5 being imaginary times, is the operator 

K =[d T - fi + h]-\ (11) 

T is the operator of kernel 

r(aj, r x \y, T y ) = K(x, r x \y, r y )K(y, r y \x, t x ), (12) 

and finally, 9(x, r) is a Gaussian field of covariance [C _1 — Y/Z\ (x, r x \y, r y ). 

We recover TF theory by taking the semiclassical limit in the Hartree energy. Indeed, in 
this limit, the density po(x) and the integrated density of states D(e) become 

Po(x) = 7d/u+ 2 (a0, (13) 

where n+(x) = (fi — W(x))§(p — W(x)), and 

D(e)= ld j dx(e-W(x)) d l 2 , (14) 

where 72 = (27r) _1 and 73 = (37r 2 ) -1 , and /i is fixed by the constraint 

Jdxp (x) = l. (15) 

At the next stage, we compute the correlation energies. For E c; i, we simply need to take 
the limit e — >• 0. For E C; 2, we first compute the terms up to the order Z~ 3 , then take the 
semiclassical limit. We recall that the true ground state energy is given by E = Z 2 ~ 2 l d E. 




,6? 0.03 
I 0.02 
0.01 

10 20 30 40 50 60 

JV 

FIG. 1: Per electron correlation energy of neutral atoms. Data points correspond to experimental 
values 



111 ] and numerical values obtained with an extended HF approach [121 ]. available for atoms 
containing up to 18 and 55 electrons, respectively. For sufficiently large atoms, the contribution 
E c /N = — 0.062 In N — 0.18 deviates from experimental and extended HF values essentially by an 
oscillating term. 

In the case of neutral atoms, and moving to standard atomic units, we find 
E c . x = -0.062 NkiN- [0.034 + 2x 2 }N + 0(N 2/3 ), 

where 

Tq being certain operator on L 2 (M. 3 ), and 

E c , 2 = 0.98 N + 0(N 2/3 ). 
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It is remarkable that the dominant term for the correlation energy —0.062 iVlniV is the 
same as the logarithmically divergent second order correlation energy of the jellium [lOj, an 



effective cut-off e on some momentum making the integrals finite in the case of atoms. Other 
parts of the correlation energy are the same as in the jellium. For its correlation energy, the 
atoms looks like a jellium with a local Fermi surface defined by p 2 < fi + (x). 

To summarize, for neutral atoms we have found the correlation contribution E c = E c -\ + 
E c -2 to the ground state energy 

E c = -0.0Q2 N\nN + cN + 0(N 2/3 ), (19) 

and not, as has been suggested, a iV 4 / 3 contribution . In figure CQ, a comparison is made 

between experimental ll|] and extended Hartree-Fock 12j values for the correlation energy, 
and formula f|T9|) taking c = —0.018 for the best fit. We have not computed numerically 



In the case of quantum dots, and in the units of ([I]), we find 

E c . x = -0.15 N + 0(N 1/2 \nN), 



(20) 



and 



E c . 2 = -0.15 N + 2 Cl N + 0(N 
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where 



ci = 
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It can be noticed that the dominant contribution to E c -i is universal, i.e. independent of 
the confining potential. 

Let us now look at the Hartree-Fock type contribution E HF = E n + E x , which needs to be 
computed semiclassically. E H depends essentially on the density of states d(e) = d e D{e). It 
is standard now to decompose d(e) into two parts, semiclassically, as d(e) = d s (e) + d osc (e) 
[l^ : and, accordingly, to decompose E UF as E HF = E HF . S + E UF - osc . The smooth part d s (e) is 
given by an asymptotic expansion in e, whose coefficients are some integrals depending on 
W(x). But the oscillating part 



is given by a sum over the periodic orbits 7 of a classical particle moving in the potential 
W(x), where S(e, 7) is the classical action along the orbit, a 7 is the orbit's Maslov index, 
and A e (e, 7) depends on the orbit's period and stability. d osc (e) therefore depends crucially 
on the nature of the dynamics associated to W(x). 

In the case of atoms, the smooth part is known up to the order e _1 , and further correc- 
tions both to the Hartree and exchange energies would be required to reach the order e of 
the correlation energy. The semiclassical determination of these corrections is particularly 
delicate due to the Coulomb singularity of W(x) at the origin, and the slow decay of the 
Coulomb potential appearing in the exchange term. Therefore, we have not undertaken this 
task. The oscillating part was computed by Schwinger and Englert, who did not give a dy- 
namical interpretation of it. The spherically symmetric potential W(x) makes the dynamics 
integrable in this case. 

In the case of dots, we have completely determined the smooth part of E KF at the same 
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FIG. 2: Oscillating contribution to the ground state energy for the chaotic quantum dot of quartic 
confinement (A = 1.2), computed with the 30 shortest periodic orbits of the classical system. 

order as the correlation energy. It is given by 

12 r N 2 f dxdy 



F, »r2 N f , / N N f dxd V / N / \ 

Ehfjs = fiN - — dxfi + (x) - — / -j -fi + (x)fj,+ (y) 

47r J Stt z J \x — y\ 



2N 3/2 
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J dxn% 2 (x) + dxAW{x)ti(ii{x))+c 2 N, (24) 



where 



dx /j 1 + 2 (x)a(x) 



c 2 = — 

/ dxfi^fifatyafa) 

a(x) being the solution of the integral equation 



(25) 



1 



a(x) = 1 



dy 



2ir J \x — y 



■0Mv))a(y); 



and where 



W(x) = V(x) + i- J 



dy 



x-y\ 



(26) 
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The oscillating corrections to the ground state energy are of order N a , and of period N 1 ^ 2 , 
where a = 3/4 in the integrable case and a = 1/2 in the fully chaotic case. Note that in both 
cases and contrary to the case of atoms, the oscillating corrections are less important than 
correlation effects. If V(x) is a harmonic well, the dynamics induced by W(x) is integrable 
and the oscillating part has been computed pjj. It is particularly interesting to consider a 
case where the dynamics induced by W(x) is fully chaotic. For this purpose we have chosen 
the quartic potential V(x) = a(xf/j3 + x\fi — \x\x\), with a = 10~ 3 and (3 = 7r/4 [ToJ, Il6 |. 
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The nature of the dynamics of W(x) depends on the couplin g co nstant A, and when A = 1.2, 
the dynamics is fully chaotic (whereas that of V(x) is not) 17] . In this case 

E nw . s = 0.14 N 2 - 0.067 N 3/2 + 0.053 N (28) 

and -Ehf ; osc5 computed by summing over the 30 shortest periodic orbits, is shown in figure[2l 
To conclude, we have found a systematic way to compute the ground state energy of 
atoms and quantum dots when the number of electrons is large. A semiclassical HF theory 
has been shown to give the dominant terms. But correlation effects have been computed 
which determine the validity of HF theory. We have also shown that shell effects result from 
the dynamics of a classical particle moving in the presence of the TF potential. This allows, 
in the case of dots, to see the impact on the oscillations of the energy of a regular or chaotic 
motion of this classical particle. 

We thank C. Plocek for useful and stimulating discussions at the early stage of this work. 
This work was supported by the Fonds National Suisse de la Recherche Scientifique. 



Electronic address: nerve. kunz@epfLch 



' Electronic address: |rico.rueedi@epfl.ch 



[1] L. H. Thomas, Proc. Camb. Phil. Soc. 23, 542 (1927). 
[2] E. Fermi, Rend. Accad. Naz. Lincei 6, 602 (1927). 

[3] J. D. Morgan, in Atomic, Molecular, & Optical Physics Handbook, edited by G. W. F. Drake 

(AIP Press, 1996), chap. 10, pp. 233-242. 
[4] E. H. Lieb and B. Simon, Adv. in Math. 23, 22 (1977). 
[5] C. L. Fefferman and L. A. Seco, Helv. Phys. Act. 71, 1 (1997). 
[6] B.-G. Englert and J. Schwinger, Phys. Rev. A 32, 36 (1985). 
[7] E. H. Lieb, J. P. Solovej, and J. Yngvason, Phys. Rev. B 51, 10646 (1995). 
[8] J. W. Negele and H. Orland, Quantum Many-Particle Systems, Frontiers in Physics (Addison- 

Wesley Publishing Company, 1988). 
[9] B.-G. Englert, Semiclassical Theory of Atoms, Lecture Notes in Physics (Springer- Verlag, 

1988). 

[10] M. Gell-Mann and K. A. Brueckner, Phys. Rev. 106, 364 (1957). 

[11] S. J. Chakravorty and E. R. Davidson, J. Phys. Chem. 100, 6167 (1996). 

8 



[12] E. Clementi and G. Corongiu, Int. J. Quantum Chem. 62, 571 (1997). 

[13] M. Brack and R. K. Bhaduri, Semiclassical Physics, Frontiers in Physics (Addison- Wesley 

Publishing Company, Inc., 1997). 
[14] C. Plocek, Ph.D. thesis, Ecole Polytechnique Federale de Lausanne (2006), li- 

brary.epfl.ch/theses/?nr=3684. 
[15] O. Bohigas, S. Tomsovic, and D. Ullmo, Phys. Rep. 223, 43 (1993). 
[16] M. Brack, S. N. Fedotkin, A. G. Magner, and M. Mehta, J. Phys. A 223, 1095 (2003). 
[17] R. Rueedi, Ph.D. thesis, Ecole Polytechnique Federale de Lausanne (2009), li- 

brary.epfl.ch/theses/?nr=4406. 



9 



